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Computer  Program  Descriptions _ 


PURPOSE: 


TROPOSPHERIC  ABSORPTION  AND  DISPERSION 
OF  MILLIMETER  AND  SUBMILLIMETER  WAVES 


TROPOSPHERIC  ABSORPTION  AND 
DISPERSION  FOR  TERRESTRIAL  OR 
EARTH-SATELLITE  PATH. 

To  compute  the  complex  retractive  index 
of  the  troposphere  for  frequencies  from  1 
to  1000  GHz:  in  particular,  the  amplitude 
and  phase  dispersion  due  to  molecular  oxy¬ 
gen  and  water  vapor  lines  and  the  attenua¬ 
tion  by  suspended  water  droplets  (radii  < 
20  junt  such  as  haze,  fog,  and  clouds. 
Extended  Basic  (HP  9845S).1 
k.  C.  Allen  and  H.  J.  Liebe  are  with  the  Na¬ 
tional  Telecommunications  and  Information 
Administration.  Institute  for  Telecommuni¬ 
cation  Sciences.  Boulder,  CO  80302. 
ASIS-NAPS  document  NAPS-03988. 

This  is  a  user-oriented  program  with  descrip¬ 
tive  comment  statements.  Inputs  are 
prompted  interactively  on  the  cathode- 
ray-tube  (CRT)  display.  Graphical  output 
is  selectively  displayed  on  the  CRT  or 
plotted  in  ink  on  paper  (e.g.  with  an  hp 
9872AJ.1  A  tabular  output  is  available 
on  the  internal  line  printer. 


LANGUAGE: 

AUTHORS: 


AVAILABILITY: 

DESCRIPTION: 


and  the  real  (dispersive)  part  follows  analogously  from 

D(/)=  2  {SF‘),+N'C  +  K-  (4) 

i 

The  summations  represent  the  contributions  of  the  resonant  line 
spectra  of  atmospheric  species  where  S  is  the  line  strength  in 
kilohertz  and  the  line  shape  factors  are  f"  and  F"  in  (GHz)"1. 
Water  vapor  lines  above  1  THz  are  lumped  into  the  continuum 
terms  A^.’  and  N'c.  The  contribution  from  suspended  water  drop¬ 
lets  are  represented  by  the  terms  N‘w  and  .V^, .  The  summations 
are  computed  for  44  02  lines  plus  29  H20  lines  below  1  THz 
(line  spectra  of  trace  gases.  03,  CO,  N20.  S02,  NHj.  etc.  are 
neglected).  Also,  weaker  H20  lines,  with  center  line  intensity 
less  than  2  X  10"  J  dB/km.  are  not  considered.  These  simplifi¬ 
cations  result  in  a  large  increase  in  computation  efficiency  com¬ 
pared  with  line-by-line  programs  such  as  the  Air  Force  Geo¬ 
physical  Laboratory  tape  )3] . 

The  Van  Vleck-Weisskopf  line  shape  factors 


/  AT 2zSHi 
v^o  /  L(yo  -/ 


-(Vq-/)5  +  T-(v0  +/)5 
)-/V+r  (v0+f)2+y2 


,  =  r  (^o  -/)  +  T&  (}>o  +/)->& _ 2  1 

LO'o-/')2 +r  (Vo+ff  +  y2  "o  J 

as  modified  by  Rosenkranz  [4]  are  used  where  p0  is  the  line 
center  frequency,  y  is  the  line  width  (both  in  GHz),  and  6  is  the 
overlap  interference  parameter.  A  line-data  base  (2]  consisting 
of  p0  and  the  spectroscopic  coefficients  a ,_5  and  bx  is 
used  with  the  meterological  conditions  to  compute  the  line 
parameters  as  follows.  For  02  in  air: 


The  three  components  for  air  of  the  complex  refractivity  N 
(in  units  of  parts  per  million) 

;V  =  A'0  F  D(/)+/Ar(/)  (I) 

namely,  the  frequency  independent  refrac..vity  N0  plus  a  refrac¬ 
tive  dispersion  D(f)  and  absorption  N"(J\  are  computed  by  this 
program.  The  frequency  /  is  in  GHz  throughout.  The  frequency- 
independent  refractivity  is  computed  using 

A'0  =  (ft  —  1  )106  =  [2.589p  +(41 .68  +  2.39jf]0  (2) 

where  n  is  the  dimensionless  refractive  index,  p  is  the  dry  air 
pressure  in  kPa.  e  is  the  water  vapor  partial  pressure  in  kPa 
(barometric  pressure  P  -  p  +  e).  and  8  =  300/T  is  the  relative 
inverse  temperature  (T  in  K).  [I  |. 

The  computation  of  D(f)  and  A'”(/)  isbased  on  the  model 
of  Liebe  [2].  The  imaginary  (absorption)  part  of  the  refractive 
index  is  computed  according  to 

»"(f)  =  S(SF"),+iV;+Jv".  (3) 


Manuscript  received  June  It,  1982;  revised  July  30,  1982.  This  work 
wjs  supported  in  part  by  the  U.S.  Army  Research  Office  under  Contract 
ARO  6-82.  See  NAPS  document  no.  03988  for  8  pages  or  supplementary 
material.  Order  from  NAPS  c/o  Microfiche  Publications,  P.  O.  Box  3513. 
Grand  Central  Station.  New  York,  NY  10163.  Remit  in  advance  in  U.S. 
funds  only  $7.75  for  photocopies  or  S4.00  for  microfiche.  Outside  the 
U.S.  and  Canada,  add  postage  of  $4.50  for  the  first  20  pages  and  $1.00  for 
each  to  pages  of  material  thereafter.  $  1 .50  for  microfiche  postage, 

1  This  is  not  an  endorsement  ot  these  items. 


S  =  a i pO3  exp  [a2(l  —  0)1 

(7a) 

y  —  o3(p  +  1.3e)00'9 

(7b) 

5  =  a4p0''s. 

(7c) 

For  H20  in  air: 

S  =  6,e03'5  exp  [b2(l  —  0)) 

(8a) 

y  =  b  j(4.80e  +  p)606 

(8b) 

6=0. 

(8c) 

The  continuum  spectrum  consists  of  A*',  the  far-wing  absorp¬ 
tion  by  strong  infrared  H20  lines,  plus  Ar‘J.  the  nonresonant  dry 
air  spectrum.  The  dispersive  part.  iYc,  which  is  nearly  zero,  is 
neglected.  Absorption  by  the  infrared  lines  is  computed  by 

A”=  l.9/ep03  1  X  10~6.  (9a) 

Equation  (9a)  is  an  empirical  expression  [S]  which,  so  far.  has 
not  found  a  generally  acceptable  theoretical  explanation.  The 
non-resonant  air  spectrum  is  computed  using 

=6.2/pe2  X  10-4  •  (l7o Hf2  +fl)l  +  2.1p8°  *  X  10-7} 

(9b) 
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where  the  lust  term  in  the  braces  is  due  to-  02  and  the  second 
to  N: ,  The  width  yQ  is  calculated  with  (7b)  using  a  spectroscopic 
coefficient  of  a}  =  0.012,  although  the  value  a3  =  0.0056  was 
found  1 1 1  using  Rosenkranz's  [4)  scheme  to  tit  laboratory 
data. 

The  liquid  water  extinction  A”  (3)  of  droplets  with  radii 
smaller  than  20  ;um  (cloud,  fog.  haze,  aerosol  hydrometeors) 
is  computed  with  the  Rayleigh  absorption  approximation 

A;;  =  4.4W/I(e  +  2)2  +  te")3j  (10) 

of  Mie  scattering  losses  where  w  is  the  liquid  water  content  in 
g/mJ.  The  )-  and  6-  dependence  of  dielectric  data  (e‘.  e") 
of  bulk  water  are  calculated  with  the  Debye  model  detailed 
in  [t>)  For  frequencies  above  300  GHz.  the  rough  approxima¬ 
tion 

a;;  =0.55u/-0  >0~6  (11) 

is  used.  The  dispersive  part  A'w  is  neglected  since  it  is  very  broad 
band  and  is  less  than  I  ps/km  for  each  I  g/mJ  of  liquid  water 
content. 

The  specific  power  attentuation  a  in  dB/km  and  the  excess 
dispersion  propagation  delay  time  t  in  ps/km  are  computed  by 

o  =  0.1S20/A  ”  (12) 

and 

t  =  3.336/).  (13) 

The  output  of  the  program  is  the  attenuation  and  excess  disper¬ 
sion  delay  for  the  chosen  path.  A  sample  output  is  given  in  Fig. 
1.  When  running  the  program,  the  minimum,  maximum,  and 
spacing  of  the  frequencies  for  which  the  refractivity  is  to  be 
computed  are  input  along  with  the  selection  of  a  terrestrial 
or  earth-satellite  path.  For  a  terrestrial  path,  the  barometric 
pressure,  relative  humidity,  atmospheric  temperature,  liquid 

water  content  in  g/m3.  and  path  length  are  input.  If  the  speci¬ 
fic  attenuation  (dB/km)  is  desired,  then  the  path  length  is  set 
to  one  kilometer.  For  an  earth-satellite  path,  the  path  eleva¬ 
tion  angle  and  ground  terminal  elevation  above  sea  level  are 
input.  An  eight-point  tropospheric  profile  of  the  meteorologic 


parameters  is  read  from  data  statements  by  subroutine  “Profilel” 
and  used  for  path  computations.  The  user  can  provide  an  alter¬ 
nate  profile.  The  cosecant  function  is  used  to  compute  the  eleva¬ 
tion  angle  dependence  of  path  attenuation  and  delay,  and  there¬ 
fore  the  program  should  only  be  used  for  path  elevation  angles 
greater  than  6°.  The  attenuation  and  delay  through  each  layer 
are  summed  to  give  the  path  total. 

Above  about  30  km.  Zeeman  splitting  of  02  lines  due  to  the 
geomagnetic  field  occurs.  Since  the  program  does  not  include 
this  effect,  output  for  earth-satellite  paths  may  be  inaccurate  for 
frequencies  closer  than  100  MHz  to  an  02  line.  The  profile  of 
meteorologic  parameters  should  be  extended  into  the  mesophere 
for  high  flying  aircraft-to-satellite  path  applications. 

Subroutine  “Gasl”  computes  the  specific  attenuation  and 
delay  from  the  meteorologic  data.  It  is  felt  that  this  routine 
along  with  the  resonance  line  data  would  be  of  interest  to  users 
desiring  to  write  programs  of  their  own  concerning  tropospheric 
dispersion. 

The  execution  time  of  the  program  requires  approximately  7  s 
for  each  frequency  and  set  of  meteorologic  conditions.  Programs 
in  Fortran  IV,  including  the  mesosphere  and  addressing  the 
Zeeman  effect  of  02  lines,  are  available  also  from  the  authors. 
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